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Abstract 

We study the interplay of the chiral and the color superconducting phase transition in an ex- 
tended Nambu-Jona-Lasinio model with a multi-quark interaction that produces the nonlinear 
chiral-diquark coupling. We observe that this nonlinear coupling adds up coherently with the o;^ 
interaction to either produce the chiral-color superconductivity coexistence phase or cancel each 
other depending on its sign. We discuss that a large coexistence region in the phase diagram is 
consistent with the quark-diquark picture for the nucleon whereas its smallness is the prerequisite 
for the applicability of the Ginzburg-Landau approach. 
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The findings of recent ultra-relativistic lieavy ion collision experiments have stimulated 
a paradigm shift, that is, quark gluon plasma (QGP) is not a weakly interacting near ideal 
gas but a strongly interactingnear perfect fluid, called sQGP, at least slightly above the 



transition temperature 



Quantum chromodynamics (QCD) exhibits a variety of 



forms of matter also at high density; chiral symmetry restoration, deconfinement, and color 



superconductivity (CSC) 



As for CSC, various sub-phases at intermediate density are 
discussed recently in addition to standard two flavor superconductivity and color flavor 
locking. 

First principle lattice QCD simulations describe high temperature phenomena but their 
applicability to finite density is limited due to the well known sign problem and/or zero eigen 
values of the fermion matrix. Therefore effective models, such as the Nambu-Jona-Lasinio 

or the random matrix model [o], [l^ , must be employed. Recently 
that handles not only chiral restoration but also deconfinement is 



(NJL) model 
the Polyakov-NJL mode 
actively studied 111, 12 



13 



Before CSC came into consideration, standard effective models predicted that chiral 
restoration at high density is a first order transition. (See also Ref. jl5] and references 
cited therein that go beyond the mean field approximation.) However, the vector interac- 
tion, which is not forbidden from symmetry consideration, rather necessary from a view 



point of nuclear physics but often ignored, may change the situation 



16 



171 ]. The competi- 



tion between the chiral (qq) and the diquark (qq) condensates on the temperature-chemical 
potential (T-/i) plane was first considered by Berges and Rajagopal 3]. In their calculation, 
the two types of condensates are mutually exclusive. But, in principle, they can coexist, 
i.e., quarks dressing chiral condensate can pair up. Pairing between such constituent quarks 



would lead to the quark-diquark picture for the nucleon 19|. On the other hand, recently 
Hatsuda et al. obtained interesting results including a new end point induced by the U(l) 
anomaly in the three flavor case, using a model-independent Ginzburg-Landau (GL) ap- 
proach to the T-/i phase diagram {20]. Here it should be noted that the GL approach is 
applicable when all the order parameters considered are small, since the free energy is ex- 
panded with respect to them. In the present case, both a and A should be small. This 
situation is realized in the vicinity of each phase transition, if the coexistence region is small 
or vanishes. It is thus an important information how large the region is, and effective models 
are useful to answer the question. 



2 



We adopt a framework that can handle these aspects of the QCD phase diagram on the 
same footing — an extended NJL model with mul ti-q uark interactions. (For introduction 
of multi-quark interactions, see also Osipov et al. 



22| and Huguet et al. |23|, [24.) In 



21 



a previous paper 



251 ]. we found that the o'^lJ^ and the interactions sharpen the chiral 



transition weakened by the uP' interaction and also that the cr^ interaction shifts the critical 
end point to a higher T, lower /i point. In the present paper, we discuss the interplay 
of the chiral and the diquark condensates brought about by the cr^A^ interaction and its 
coherence/competition with the iJ^ interaction. Here, cr, tu, and A denote the scalar, vector, 
and diquark auxiliary mean fields, respectively, defined later. Note that our study is limited 
to the two flavor case at present. 

The Lagrangian density of the extended NJL model adopted in the present work is given 

by 



q{% ^ - mo)q+ ^g2fi {^{qq) + {qi'^->Tq) j +5(4,0 \{qq) + {qil-,Tqy 

-goMi'^qf - 92,2 ((gg)' + (g^Tsrg)') (57''?)' 

+do,2({ifee^-f5q){iqee^-f5q^) + {iq^ee''-f5q){iqee''l5ql 



+d2,2{qqf[{iq^ee^-i^q){iqee^-i^q^) + {iq^ee^-i^q){iqee^-i^q^) ) + 



(1) 

where q'^ = Cq^ and (f = q^C are the charge-conjugation spinors, C = i'-f'^'-f^ is the charge- 
conjugation matrix, q is the two flavor quark field, r = (r^,r^,r^) are the Pauli matrices, 
mo = diag(m„,m(i) is the current quark mass matrix, e and are the totally antisymmet- 
ric tensors in the flavor and color spaces, and Qij and d^^n {hji rn,n = 0, 1, 2, ■ ■ ■ ) are the 
coupling constants of quark-quark interactions. We consider only the four- and eight-quark 
interactions ignoring higher-order interactions denoted by ellipsis in the Lagrangian den- 



sity. Among the physically important eight-quark interactions, 54,0 and 52,2 terms wil 



ignored in the following, since their roles have been clarified in our previous paper 25| as 



be 



mentioned above and we would like to concentrate on the chiral-diquark coexistence. Under 
the standard mean field approximation (MFA), the Lagrangian density reads 

Cmfa = g(^ ^ - (mo + S3) + Kl')q - ^^ti^l'^^'l^q) - ^^'^{zqee'-f,q') - U, (2) 
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where 



yb 



-2{g2,Qa + (i2,2cr|A|^), 
-25^0,2^^, 

-2(4,2A*' + rf2,2^'A*^), 

U = ^2,0^7^ - 50,2^' + cio,2A*''A'' + 3ci2,2a2 A*''A^ (3) 

and the auxihary fields introduced are a = {qq), to = (g7°g), A^ = {iifet^^^q), and A**" = 
{iq^ee^-ir^q). 

The thermodynamical potential Vl of the system with finite temperature T and chemical 
potential is then obtained as 



+sgn{E;)E^ + E+ + |{ln(l + e-^g°(^p )/^^a) + ln(l + e-^^^)}] 
+VU, (4) 
where /3 = l/T, /i = // + E^, M = mo + S^, = ^p^ + M^, ^ ± jl, E^ = 



-E^^ + lEdp, and sgn(£'p) is the sign function. The corresponding scalar, vector, and 
scalar diquark densities, Ps, py and p^ are given by 

Ps = -2NfM I ^-l{l-n(£;p-)-n(£;+) + ^tanh^ + Stanh^}, (5) 



■ ''^'tanh^}. 
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2 £;+ 2 J' 



where 



no 



(6) 
(7) 

(8) 



" l + exp{/3(£;p-/i)}' '"^ l + exp{/3(£;p + /i)}' 
The gap equation can be derived by minimizing the thermodynamical potential with respect 
to cr, a;, and A*, their physical solutions then satisfy the stationary condition 

/ d_ (n\ \ 
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where the effective couphngs are 

/ 
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c c c 
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aA* 




aA* 
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2(^72,0 


+ C^2,2|A 


^) 






V 





4d2,20-A 



25-0,2 




(10) 



4c?2,2CtA* 


2(rfo,2 + d2,2CT^) 

When det(^) ^ 0, Q has its inverse, and then the stationary condition leads to cr = ps, 
oj = Pv, and A = p^. 

It has been shown that the effect of the couphng on the phase diagram is suppressed by 
the non-hnear terms, QAflO* and (72,20"^"-^^; in our previous paper j25|. Therefore, the vector 
couphng 50,2 is fixed to the small value, 0.252,0- The adopted parameters for numerical 
calculations are summarized in Table I. We examine both signs for ^2,2 since this is not 
determined within the model and they would lead to different physical pictures. 



model 


92,0 


50,2 do,2 




NJL + A2 


5.498 


Ga 





NJL + + A2 


5.498 


Gui Ga 





NJL + A2 + (j2A2 


5.498 


Ga 


0.2Ga 


NJL + a;2+A2 + a^A^ 


5.498 


Guj Ga 


0.2Ga 



TABLE L Summary of the parameter sets. The coupling constants are shown in GeV ^. For all 
cases we adopt mo = 0.0055 GeV, A = 0.6315 GeV, and o^^ = -0.03023 GeV. Here, G^ and Ga 
are 0. 252,0 and 0.6^2,0) respectively. (See Table 1 in Ref. [25] for comparison.) 



In the following, we discuss the phase diagrams obtained by adopting the models with 
the parameters summarized in Table I putting emphasis on the chiral-diquark coexistence at 
low-T. Figure [U^a) graphs the phase diagram of the standard NJL model with the diquark 
condensate. The coexistence region is very small in this case. Blaschke et al. first pointed 
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out the existence of this coexistence region adopting another parameter set without the 



lJ^ interaction 26|] . Figure [11(b) shows the effect of the <J^ interaction, that is, it weakens 
both transitions and shifts the chiral restoration to the higher density side and consequently 
produces a coexistence phase at low-T. This confirms the result presented by Kitazawa et 
al. 171]. Comparison of Figs. [11(a) and [21(a) demonstrates the effect of the a^A^ coupling for 
the case of positive (i2,2- This nonlinear interaction shifts the CSC transition to lower density 
and consequently produces a coexistence phase. Figure [2]|^b) includes both the uP' and the 
cr^A^ interactions. They coherently add up in this case. This result can be understood from 
the expressions of the effective couplings that lead to 

Ss = -Gs.Or=-2((72,0 + C^2,2|A|2)or (H) 

and 

Ed = -GdA A = -2 (rfo,2 + ^2,2^') A, (12) 

in which the former indicates that positive ^2,2 enhances |Ss| when A 7^ exists and the 
latter indicates that positive c/2,2 enhances |Sd| when a 7^ exists. This means that the a^A^ 
interaction acts only when the lo^ interaction makes a and A coexisting. Figures [3]^a) and 
[ni^b) graph the result of the negative ^2,2- In this case the lJ^ interaction and the a^A^ one 
are destructive to each other and consequently the coexistence region becomes very small. 

As mentioned above, the sign and the magnitude of the eight-quark interaction is not 
determined within the present model. The positive ^2,2 results in a large coexistence region. 



The sign is supported by the result of the quark-diquark model for the nucleon [19| that 
the diquark interaction is sizably stronger in the normal baryon-number density region than 
in the high density one. However, further analysis is needed to determine the strength of 
the coupling more precisely in the normal density region. Oppositely, if the negative sign is 
favored by some reason, the coexistence region shrinks and the two phase transitions occur 
at almost the same T and /i. In the vicinity, the order parameters, a and A, are small. Such 
a region is an ideal playground for the GL approach, since the free energy is expanded with 
respect to them. In other words, when the sign is positive, the GL model is not useful to 
determine the phase diagram except for the high-T region where both the order parameters 
are small. 

We have studied the interplay of the chiral and the color superconducting phase transition 
in an extended Nambu-Jona-Lasinio model with a multi-quark interaction that produces the 
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FIG. 1: Phase diagrams given by (a) the standard hnear NJL model and (b) the extended NJL 
model that includes the interaction. 
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FIG. 2: Phase diagrams given by (a) the extended NJL model that includes the cj^A^ interaction 
and (b) that includes the and the cr^A^ interactions. In these calculations ^2,2 is positive. 

nonlinear cr^ coupling. We have found that the size of the chiral-diquark coexistence region 
is sensitive to the sign of the coupling. The positive sign is supported by the quark-diquarlc 
model for the nucleon, but further analysis is needed to determine the density dependence 
of the diquark interaction more precisely. Meanwhile, the negative sign is the prerequisite 
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FIG. 3: The same as Fig. [2] but ^2 2 is negative. 



for the appHcabihty of the Ginzburg-Landau approach that has already been apphed to 
determine the phase diagram. Thus, the determination of the sign is an important subject 
related to the phase diagram. 
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